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Layered doped Mott insulators, such as the cuprates, show unusual temperature dependence of 
the resistivity. Intriguingly, the resistivity perpendicular to the CuC>2 planes, p c (T), shows both 
metallic (dp c /dT > 0) and semi-conducting (dp c /dT < 0) behavior. We shed light on this puzzle by 
calculating p c for the two-dimensional Hubbard model within plaquette cellular dynamical mean- 
field theory and strong-coupling continuous-time quantum Monte Carlo as the impurity solver. 
The temperature, T, and doping, S, dependencies of p c are controlled by the first-order transition 
between pseudogap and correlated metal phases from which superconductivity can emerge. On the 
large doping side of the transition p c (T) is metallic, while on the low-doping side p c (T) changes 
from metallic to semi-conducting behavior with decreasing T. As a function of doping, the jump 
in p c across the first-order transition evolves into a sharp crossover at higher temperatures. This 
crossover coincides with the pseudogap temperature T* in the single-particle density of states, the 
spin susceptibility and other observables. Such coincidence in crossovers is expected along the 
continuation of the first-order transition into the super-critical regime, called the Widom line. This 
implies that not only the dynamic and the thermodynamic properties but also the DC transport in 
the normal state are governed by the hidden first-order transition. p c (T) has a high-temperature 
quasi-linear regime where it can exceed the Mott-Ioffe-Regel limit and when it has a minimum it is 
nearly parallel to the Widom line. 

PACS numbers: 71.27.+a, 71.10.Fd, 74.25.-q, 71.30.+h 



The puzzling behavior of electrical resistivity has been 
at the center of the high-temperature superconductivity 
conundrum from the very beginning ' . Normal-state re- 
sistivity is highly anisotropic : the out-of-plane c-axis 
resistivity p c {T) can be orders of magnitude larger than 
the in-plane resistivity p a b(T). In addition, for values 
of doping where the cuprates become superconducting, 
the normal-state p a b(T) and p c (T) can show contrast- 
ing behaviors: the in-plane resistivity p a b(T) is metallic, 
while the out-of-plane c-axis resistivity p c {T) can be both 
metallic or non-metallic (semi-conducting) . The opening 
of the pseudogap can be detected by the deviation in 
both the in-plane and the c-axis resistivities from their 
linear-T behavior at high temper atures 5-7 . 

The behavior of the c-axis resistivity, its relation to 
superconductivity and to the pseudogap are important 
issues for the understanding of cuprates. Diverse expla- 
nations for the interplane resistivity have been offered 
over the years. They can be classified into two groups, 
depending on whether it is unconventional for the in- 
plane physics or for the interlayer coupling. The first 
group includes mechanisms based on spin-charge separa- 
tion 8 ' , fluctuations of the phase of the superconducting 
order parameter 10 or in-plane strong-coupling physics . 
The second group contains models where the interlayer 
tunnelling is influenced by disorder 12 , by bosons 13-15 or 
by interplane and in-plane charge fluctuations 16 . 

Here we examine the interplay between c-axis trans- 
port, superconductivity and pseudogap by studying the 



DC c-axis resistivity of a hole-doped Mott insulator rep- 
resented by the one-band Hubbard model. We explore 
the possibility that Mott physics -essentially the blocking 
of charge motion driven by strong Coulomb repulsion- 
can determine the intricate doping and temperature be- 
havior of p c . In our approach, the interplane transport is 
governed by the in-plane scattering, as in the first group 
of mechanisms for unconventional interplane conduction. 

We solve the model using cellular dynamical mean- 
field theory (CDMFT) 17 ' 1S for a self-consistent 2x2 pla- 
quette. Recent developments in the algorithms make 
the present study possible. Recent work 11 using a 2-site 
cluster strengthen the experimental correlation between 
the behavior of p c (T) and the opening of the pseudogap 
as revealed by angle-resolved photoemission (where the 
pseudogap appears as a lack of a quasiparticle peak at the 
antinode ) and c-axis optical conductivity cr c (w) (where 
the pseudogap appears as a low-frequency suppression of 
spectral weight transferred to high frequencies 21 ). Sim- 
ilar c-axis optical conductivity results have been shown 
in Ref. 22 using an 8-site cluster. 

A recent development requires that those findings be 
re-examined: in plaquette CDMFT, a first-order tran- 
sition ' occurs at finite doping between a pseudogap 
and a correlated metal. This transition is connected to 
the Mott transition in the undoped model. The crossover 
to the pseudogap state, T*(S), lies along the thermody- 
namic crossover (known as Widom line ) that begins at 
the critical endpoint of the transition and extends in the 
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supercritical region. This indicates the common origin of 
the pseudogap and of the thermodynamic crossovers 21 '. 
The further step provided by the present work is to ad- 
dress the role of the Widom line for interplane transport. 
In agreement with previous studies '", we attribute the 
semi-conducting p c {T) to the development of the pseudo- 
gap. We further show that the temperature and doping 
dependence of p c is governed by the Widom line crossover 
generated by the pseudogap to correlated metal first- 
order transition, thereby providing a unified picture for 
explaining DC transport, thermodynamic and dynamic 
properties. The resistivity minimum, Ioffe-Rcgcl limit, 
inflection points and their relation to the Widom line are 
also discussed. 

Method- We consider the two dimensional Hubbard 
model on a square lattice, 



vertex correction and the c-axis conductivity is given by 
the bubble diagram 
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with tij the nearest neighbor hopping amplitude, /x the 
chemical potential and U the on-site Coulomb repulsion. 
The operators c^. and Ci„ respectively create and annihi- 
late an electron with spin a at site i and i%i C = cJ.Ci CT 
is the number operator. We solve this model using 
CDMFT 1 '' 18 ' 2 ' . In this approach, a cluster of lattice 
sites, here a 2 x 2 plaquette, is embedded in a self- 
consistent bath of non-interacting electrons. The action 
of the cluster coupled to the bath is given by 
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S = S C + dr dr'^(r)A(T,T')V(T'), (2) 



with S c the action of the cluster and A the hybridization 
matrix. We can then self-consistently determine A by re- 
quiring that the infinite lattice and the plaquette have the 
same self-energy and the same Green's function on the 
plaquette: A(iuj n ) = iuj n + p — t c — E c (iw„) — G{i(jj n )~ l . 
Here t c is the cluster hopping, E c is the cluster self- 
energy, and G(iu) n ) = J^k 
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superlattice wave vector. To solve the impurity prob- 
lem Eq. (2), we use a continuous-time quantum Monte 
Carlo method 19 ' 28 ' 29 which sums all diagrams obtained 
by the expansion of the action Eq. (2) in the hybridiza- 
tion A. With this approach, we constructed the normal 
state phase diagram of the model in Refs. 23,24,26, and 
the superconducting phase diagram in Ref. 30. 

Here we consider the c-axis DC electrical conductiv- 
ity. We work in the limit where the hopping t± along 
the c-axis, i.e. between planes, is much smaller than in- 
plane energy scales, t and U. To second order in t± 
then, the Green's functions entering the Kubo formula 
for the conductivity can be evaluated to zero'th order in 
t± since each one of the vertices is already of order t±. 
In addition, since t± at the vertices create an electron- 
hole pair where the electron and the hole lie on different 
planes, they do not interact since the Coulomb repulsion 
U is local. Therefore, to leading order in t± there is no 
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The factor 2 in front comes from the spin summation, 
(7° = e 2 c/hab, with a, b in-plane and c interplane lat- 
tice constants, f(u>) is the Fermi function, A(k, ui) — 
— — ImG(k, u;) is the one-particle spectral function, and 
t± (k) = £o(cos k x — cos k y ) 2 is the interplane hopping for 
collincar Cu02 planes' 52-34 . Clearly, the k dependence of 
t_i_(k) makes the c-axis conductivity mainly sensitive to 
the antinodes, where the pseudogap develops . 

To compute the DC conductivity, one has to 
take the limit SI -) 0. We further simplify 
/(a,)-/(*+n) ^ _<MM ^ g^) to obtain a c {Vt = 0) = 

2cr° J^k t (k)^! 2 (k, u = 0). This expression is convenient 
because a simple quadratic extrapolation of the Matsub- 
ara data for the Green's function towards lo = suffices 
to obtain the spectral weight at the Fermi energy. This 
expression for the conductivity is a good approximation 
if the spectral function does not vary much for |cj| < T. 
We have verified that this is a good approximation. The 
last step is to construct the lattice A(k,oj = 0) from 
cluster quantities. There are different methods to do this 
interpolation. Here we use the Green's function peri- 
odization , which follows directly from the Self-Energy 
Functional Approach 2 ' ' 5 ' . We checked that cumulant pe- 
riodization 58 gives qualitatively similar results. 

We consider the following convenient model parame- 
ters: U — 6.2t, which is larger than the critical value 
of the Mott transition at half-filling, Umit ~ 5.95 23,39 , 
t = 1, to = 0.05i and we present p c data in units of 
p° = 1/(7°. At larger values of U, the critical point of the 
first-order transition moves to dopings that are larger 
and more consistent with cuprate phenomenology, but it 
is then located at temperatures too low to be accessi- 
ble numerically 21 . To convert into physical units we use 
t = 0.35eV and c/ab = 0.5 A" 1 . 

Temperature dependence of c-axis resistivity - Figure 1 
shows the normal state c-axis resistivity p c as a func- 
tion of temperature T for different values of hole dop- 
ing 5. The insets show the computed superconducting 
mean- field temperature T c d , below which Cooper pairs 
form within the 2x2 plaquette ! "' J " if we do not allow 
for antiferromagnetism. Buried below the superconduct- 
ing phase (blue/light grey region), there is the first-order 
transition from a pseudogap state to a correlated metal 
terminating in a critical point at finite T and finite 8 
(red/dark grey region). 

We study p c (T) in a broad range of temperature: at 
high T, p c {T) has an approximate linear behavior with 
no sign of saturation. It can exceed the generalized Mott- 
Ioffe-Regel maximum metallic resistivity saturation limit 
for anisotropic systems, p a bPc < p\- This has also been 
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FIG. 1: (Color online) c-axis resistivity p c /Pc versus temper- 
ature T for several values of hole doping 5 = 1 — n. Data are 
obtained within CDMFT on a 2 x 2 plaquette in the normal 
state for U = 6.2. Thick segments mark the superconducting 
transition temperature T^. Inset: (S, T) phase diagram. A 
first-order transition (red/dark grey region) between a pseu- 
dogap and a correlated metal terminates at a critical point 
(red dot). It lies below the superconducting phase (blue/light 
grey area). Vertical lines indicate the values of doping of the 
resistivity data in the main panels, a) For dopings above 
the first-order transition, p c (T) decreases monotonically with 
T. b) For dopings below the first-order transition, p c (T) has 
a non-monotonic behavior. Filled dots mark the resistivity 
minimum. 



observed in DMFT for isotropic systems, as discussed 
recently 41 " 43 . At low T, CDMFT allows us to study the 
underlying normal state of the model by suppressing the 
superconducting order. Hence, we can reveal the normal 
state p c below T%. 

The overall magnitude of p c decreases with increas- 
ing doping and presents two characteristic behaviors. (1) 
For values of doping larger than the first-order transi- 
tion (Fig. la), p c increases monotonically with T, i.e. 
has a metallic behavior, dp c /dT > 0. Moving away 
from the first-order transition, the zero-temperature in- 
tercept p c (T 0) decreases, and in the low-T regime 
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FIG. 2: (Color online) Normal-state c-axis conductivity a c /cr° 
versus doping. For T <T P , where T p is the temperature of the 
second-order critical endpoint, a c (8) has a discontinuous jump 
(line with circles). Above T p , a c (S) is continuous and rapidly 
crosses over from a low interplane conductive state (the pseu- 
dogap, PG) to a high conductive state (the correlated metal, 
CM). The inflection point in a c (p) defines the characteristic 
crossover temperature T„ c (S). Inset: —d(a c /a®)/dp versus 5 
in a semilogarithmic scale. The maximum of the peak is used 
to locate the inflection point in a c (p). 



p c {T) gradually approaches Fermi liquid behavior (see 
<5 = 0.25). (2) For values of doping smaller than the first- 
order transition (Fig. lb), p c (T) has a non-monotonic 
behavior. p c (T) is metallic at high temperatures, goes 
through a minimum at T m ; n (full circles) and crosses over 
to a semi-conducting T dependence (dp c jdT < 0) as T is 
decreased. T min decreases with increasing doping, lead- 
ing to a wider range of T with metallic dependence when 
5 is increased toward the first-order transition. 

Why does p c (T) increase with decreasing T? Many 
works have rooted the upturn in p c {T) to the onset of 
the pseudogap state 1 ,J 1 . Within cluster DMFT methods, 
previous work reported both the occurrence of a pseudo- 
gap state close to the Mott insulator 26 ' 27 ' 35 ' 44 , and the 
correlation between the semi-conducting-like p c (T) be- 
havior and the pseudogap phase . Our systematic data 
in Fig. 1 confirm these findings, but recast these results 
in a new framework: the first-order transition between 
the pseudogap and the correlated metal defines a wa- 
tershed for transport properties. It separates a regime 
where p c (T) shows a metallic behavior from a regime 
where p c {T) has a non-monotonic behavior. 

Doping dependence of c-axis resistivity.- To further 
characterize how the first-order transition modifies the 
c-axis transport, we present in Fig. 2 the c-axis conduc- 
tivity cr c as a function of doping for several temperatures. 
For T <T p , a c (S) is discontinuous at the transition with 
a jump of almost one decade of magnitude (lines with cir- 
cles). Coexistence is found for a c (p) in the (T,p) plane. 
The jump in the conductivity is a hallmark of the first- 
order nature of the transition. 
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FIG. 3: (Color online) Temperature versus doping phase dia- 
gram of the two dimensional Hubbard model within plaque- 
tte CDMFT for U = 6.2. Below the superconducting region 
delineated by (blue/light grey area), the first-order transi- 
tion (red/dark grey area) terminating at the critical endpoint 
(<5 P , T p ) (circle) separates a correlated metal from a pseudogap 
metal. T aa (S) is the temperature where o" c (m) has an inflection 
point. It follows T* and Twl, i.e. the dynamic and thermo- 
dynamic supercritical crossovers determined by the inflection 
in the local density of states A(u) = 0, T) and in the charge 
compressibility n(fi) respectively. The pseudogap scale can be 
identified also as inflection points in the local spin susceptibil- 
ity Xo(T), T xo . r Pc>m i n is the temperature where p c (T) has a 
minimum. It scales with the temperature where A(u> = 0, T) 
[Xo(r)] peaks, lX ma x [r xo ,max], and can be used as a predic- 
tor of the crossover T„ r . 



At the critical temperature T p , the two distinct phases 
separated by the transition (pseudogap and correlated 
metal) merge into one. o~ c {8) is continuous with an in- 
finite slope at S p , do~ c /dp\$ p — > oo. The endpoint of the 
first-order transition generates a characteristic crossover 
in the supercritical region T > T p : the divergence in 
do~ c /dp, is replaced by a peak, whose maximum decreases 
away from T < T p (see inset of Fig. 2). From the con- 
ductivity data, we define a characteristic temperature T CTc 
as the point at which the derivative da c /dp reaches its 
maximum. T a - a (8) characterizes the pseudogap to corre- 
lated metal crossover. In Fig. 3 we show the crossover 
line T CTc (<5) (line with full orange triangles) that sharpens 
towards the first-order transition. This crossover sticks 
out of the superconducting region and can be used to 
predict the first-order transition hidden by superconduc- 
tivity. The rapid rise of conductivity (or, equivalcntly, 
drop of resistivity) with doping above the superconduct- 
ing region is a stringent prediction of our theory. 

Signature of the Widom line in the c-axis transport - In 
a recent work 2 '' we have linked the dynamic crossover cor- 
responding to the opening of the pseudogap in the density 
of states to the thermodynamic continuation of the first- 
order transition in the supercritical region, called Widom 
line. Theoretically, the latter is defined as the line where 



the maxima of different thermodynamic response func- 
tions merge close to the critical endpoint . The Widom- 
line crossover governs DC transport properties as well. 
This is shown Fig. 3 that compares three crossover lines: 
T (Jc , i.e. the interplane transport crossover in Fig. 2, T*, 
i.e. the dynamic crossover signaling the pseudogap phase, 
obtained from the inflection point in the local density of 
states 2 ' 1 A(uj = 0) along paths at constant <5, and Twl, 
i.e. the thermodynamic crossover identified by a peak 
in the charge compressibility k = l/n 2 dn/ dp (or by the 
inflection point in the T dependent spin susceptibility 
(Knight shift), T Xo ). The three phenomena are concomi- 
tant. All crossover temperatures decrease with increasing 
doping and they end at the critical endpoint (5 p ,T p ). 

Another important feature of p c (T) is the appearance 
of minima in Fig. fb. In Fig. 3 we also show the loca- 
tion of these minima Tp cjmm . While T Pcjm i n occurs at 
higher T than T CTc , the lines are nearly parallel. Thus 
7p ci mi n can be used as a high temperature predictor of 
the crossover T tJc (or, equivalently, T*). It also scales 
with Ta .max, the maximum of the local density of states 
A(uj = 0,T). However, T Pcim - m does not end at the crit- 
ical point (5p,T p ); the minimum of p c (T) becomes more 
shallow with increasing 6 (see Fig. lb) and eventually dis- 
appears for values of doping larger than the first-order 
transition. Therefore our systematic analysis rules out 
interpretations based on a linear extrapolation of T Pcim i n 
(or TA, m in) to T 0. This would lead to a value of crit- 
ical doping S R3 0.07, in contradiction with the metallic 
like p c (T) above 8 « 0.05. 

T a - c (S) intersects the superconducting phase delimited 
by T£ . This result supports our discovery that within 
cluster DMFT, the pseudogap and superconductivity are 
distinct phenomena , a result confirmed by larger clus- 
ter calculations 40 ' 45 . Here, we find that superconductiv- 
ity can appear from a normal-state where out-of-plane 
conduction is semi-conducting like (i.e. non-metallic on 
the small 5 side of the transition) or metallic (on the large 
6 side of the transition). 

In summary, the first-order transition ending at a crit- 
ical point (S p ,T p ) and its associated crossover in the su- 
percritical region emerges as the unifying mechanism to 
interpret the out-of-plane transport. It is a watershed 
separating a regime where p c (T) shows a metallic be- 
havior from a regime where p c {T) has a non-monotonic 
behavior. The rapid increase of a c with doping coincides 
with the pseudogap temperature T* and with the Widom 
line, namely the thermodynamic crossover generated by 
the first-order transition in the supercritical region. The 
resistivity minimum is distinct from T* but follows a line 
that is a large T precursor. Thus we ascribe DC trans- 
port, dynamic and thermodynamic crossovers to a com- 
mon origin. 
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